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ABSTRACT 

We show that, in the presence of the steep temperature gradients 
characteristic of EUV models of the solar transition region, the electron 
and proton velocity, distribution functions are non-Maxwellian and are 
characterized by high energy tails. We estimate the magnitude of these 
tails for a model of the transition region and compute the heat flux to be 
a maximum of 30% greater than predicted by colli si on- dominated theory. 
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1. INTRODUCTION 

Existing models of the solar transition region, derived from 
observations of the solar spectrum in the extreme ultraviolet (EUV), 
make use of the assumption that the electron velocity distribution 
functions are nearly Maxwellian to compute excitation, ionization, . and 
recombination rates. However, these models are characterized by steep 
temperature gradients which, as we shall show, lead to distribution 
functions for electrons and protons which are non-Maxwell i an in contra- 
diction with the initial assumption. 

Departures from the Maxwellian distribution can be substantial 
(depending on the steepness of the temperature gradient) and take the 
form of high energy tails for those particles which are propagating down 
the temperature gradient. This leads to heat fluxes differing from 
those computed from collision-dominated theory and ionization, excitation 
and recombination rates differing from those computed assuming Maxwellian 
distribution functions for the electrons. 

In this paper, we discuss the kinetic effects which give rise to 
high energy tails in the presence of steep gradients and estimate the 
magnitude of these effects for Dupree’s (1972) EUV model of the transition 
region. We then compute various moments of the electron distribution 
functions and compare these to results obtained from collision dominated 
theory. A subsequent paper will compare the ionization equilibrium 
populations for ions of carbon, nitrogen, and oxygen, computed with the 
non-Maxwellian distribution functions obtained in this paper to those 
computed with the usual assumption of a Maxwellian distribution function 
for the electrons. 
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2. NON-MAXWELLIAN ELECTRON AND PROTON DISTRIBUTION FUNCTIONS 
A. Departures from Maxwellian 

The collision dominated kinetic theory of non-uniform gases was 
developed by Chapman and Enskog in the 1920' s and is expounded by Chapman 
and Cowling in "The Mathematical Theory of Non-Uniform Gases" (1970). 

The Chapman-Enskog theory assumes that the distribution function for a 
given species, a, can be written in general as, 

f a (x,v,t) = f a (0) (x,v,t} + f a (1) (x,v,t) + f a (2) (x,v,t) + . . . (1) 


where each of the f ^ represents a successive approximation to the 

ct 

total distribution function and each satisfies the condition. 


f (n) 

ct 

f (n-1) 

ct 


« 1, for all x, v_» and t. The full Boltzmann equation 


is then solved with what amounts to a perturbation scheme. The zeroth 
order solution, f is a Maxwellian distribution function characterized 
by a temperature, T, density, n, and mean flow, v. 

o 

Spitzer and Harm (1953) computed the first order electron distribution 
( 1 ) 

function, f ' 1 , in the presence of a temperature gradient for a gas of 
6 

mean ionic charge 2 defined by. 


Z 1 I n.Z?/n (2) 

.lie 

with the sum taken over all positive ions. In Spitzer and Harm's notation, 
the total electron distribution function, to first order,' is given by. 
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( 3 ) 


f e = f e <0) (l + vD e Cu,z)]. 

where u is the electron speed normalized to the thermal speed and y 
is the cosine of the polar angle in velocity space. 

In Table II of Spitzer and Harm the quantity Z D g (u,z)/B e is 

[ 2k 2 T (vT j . 

tabulated as a function of u, where B j = 1 is the ratio of 

6 \ Trn e e 4 £nA 

the electron mean free path to temperature scale length and A is the 
number of particles in a Debye sphere. We note that D(u,z)/B is a 
function of u only and that this quantity increases rapidly with u. In 
Figure 1 5 we plot B as a function of log T for Dupree's transition- 
region model. 

Since departures from Maxwellian are large when |yD e (u,z)| >_ 1, we 

2 

define a critical energy (u ) such that |D (u ,z)| =1. In Figure 2 

v* c L 

2 

we plot u c as a function of log T throughout the transition region for 
Dupree's model. We find that, where the temperature gradient is steepest 
(T 'v ID c K), substantial deviations from Maxwellian exist at energies 
greater than or equal to approximately 6 kT. We will show that this 
result has important implications for ionization equilibrium and heat 
fluxcalculations. 

A similar analysis can be carried out for the protons -.given the first 
order distribution function computed by Roussel-Dupre (1979). In 
Figure 3, we plot the critical energy (u ) at which |D. (u ,z)| = 1 as a 

w I 

function of proton temperature, assumed to equal the electron temperature. 

C O 

We find that, where the temperature gradient is steepest (T ^ 10 K), 

substantial deviations from Maxwellian exist again at energies greater 
than or approximately equal to 6 kT. We now ask' in what way we may 
expect these distribution functions to differ from a Maxwellian. 
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B. High Energy Tail 

Consider a plane-parallel layer of gas, with z as the vertical 
coordinate, composed of field particles having a Maxwellian distribution 
locally with temperature, T^, and density, n^, each varying as a function 
of z according to Dupree's model for the transition region. We now 
inject test particles with velocity, v_, parallel to the z-direction and 
determine the rate at which the test particle's velocity and kinetic 
energy change as a result of encounters with field particles. From these 
rates, we then obtain the time, t^, for a test particle to undergo a 90° 
deflection and v the time, tp, for a test particle to thermal ize to the 
local field particle temperature. This problem was worked out by 
Spitzer (1962), cf., his page 132. 

Given that the field particles have a non-uniform temperature in 
the z-directi'on, we can define another time scale; namely, the time it 
takes for a test particle with velocity, _v, to travel a temperature 
scale length. We have. 


t 


F 




dz J 


(4) 


• Examination of the energy dependence of the ratios t Q /tp and tp/tp 

reveals that both ratios increase rapidly with test particle energy. It 

is possible then to define critical test particle energies, (u c )^ and 
2 

(u ) E , above which tp/tp and tp/tp, respectively, become greater than one. 

2 

In Figure 4, the solid and dashed curves represent plots of (u c )p and 
2 

(u )p, as a function of Log T through the transition region with 
electrons as the test particles. We find that, where the temperature 
gradient is steepest, the deflection time becomes on the same order as a 
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time of flight (tp) for energies greater than or approximately equal to 10 
kT. Since test electrons with energies greater than (u c can penetrate 
a temperature 'Scale length without undergoing a 90° deflection, high 
energy electrons from the corona can penetrate quite far into the 
transition region while the colder test electrons from the chromosphere 
will not penetrate as. far. This leads to an anisotropic velocity distri- 
bution function. 

Similarly, we find that t £ /tp 1 for energies > 5 kT at T = 10 5 ’^ 
and since test electrons with energy greater than (u c will maintain 
their energy over scales greater than a temperature scale length, high • 
energy electrons from, the corona will populate the high energy tail of 
the field electron distribution function. Thus, the final self-consistent 
distribution functions will possess high energy tails which are highly 
anisotropic. 

2 2 

The values of (u £ )p and (u c for proton test particles are plotted 

as functions of log T in Figure 5. The protons clearly will also form 

non-Maxwell ian distribution functions, just as the electrons. An interesting 

difference, however, stems from the fact that for the most part electrons 

and protons do. not exchange energy. A proton test particle may thermalize 

to a temperature which is completely different from the electron temperature. 

m 

This follows since t £ (p-e collisions) ^ ^ t E (p-p collisions). 

The ions responsible for the EUV line emission present a different 

scenario. Since the ion abundances relative to the protons and electrons 

are quite small, and since their masses are much larger than the electron 

mass, they will interact primarily with protons. Furthermore, the 

5 6 

average charge of ions formed' from 10 K to 10 K ranges from Z'= 3 to Z = 

2 

10 and since t Q and t £ are both inversely proportional to Z , the ions 
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will usually isotropize to nearly Maxwellian distribution functions and 
thermalize to the proton temperature. To illustrate this result, we 
have carried out an analysis on Si IV similar to that carried out for 
the electrons. In Figure 6, we plot {u c )p and (u c )^. vs. Log T, for 
Si IV and find that these critical energies are much larger than for the . 
electrons and protons. Thus the ion distribution functions are nearly 
Maxwellian. An interesting consequence of these results is that EUV 
line profiles will reflect the proton temperature but not necessarily 
the electron temperature and will not reflect any asymmetries which 
might arise in' either the electron or proton distribution function. We 
note, however, that the latter conclusion applies only to the extent 
that we have neglected the first order correction (to the zeroth order 
Maxwellian distribution function) obtained from collision-dominated 
theory. - Indeed, thermal diffusion (a process which is manifested in the 
first order ion distribution functions [cf. , Roussel-Dupre,1979] ) 
of the heavy ions clearly reflects the asymnetries associated with the 
first order electron and proton distribution functions. 

As a final note we emphasize that, while it is true that there will 
always exist particles with high enough energy to be collisionless in 
the transition region, the supply of such particles is limited by the 
ultimate source of energy, i.e., the solar corona. Thus, let us assume 
that the coronal electron distribution function is a Maxwellian with 
temperature 10 K, the fraction f of particles with energy greater than 
u 2 is listed in Table I. We find that at T = 10 5 K, for example, 76 
percent of coronal electrons have energy greater than 6 kT; however, 

.only two percent have energy greater than 50 kT. The point is that the 
critical energies at which electrons become collisionless are sufficiently 
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low through the upper half of the transition region that a substantial 
number of coronal electrons are collisionless through that region. As a 
result,, large departures from Maxwellian' can be expected in the region 
from 10 5 to 10^K. 


3. PARTICLE ORBITS ~ 

In the previous sections we showed that high speed test electrons, 
moving along the temperature gradient in the transition region, could v 
penetrate a temperature scale length without undergoing a 90° deflection 
or losing a significant fraction of their initial energies. The analysis 
which brought us to these conclusions, however, was based only on a 
rather simple minded comparison of appropriate time scales. .In this 
section, we shall develop these ideas on a more quantitative basis. 

We shall first .discuss some properties of the ’dynamical -friction 1 
force exerted on a test electron as it moves through a fully ionized, 
gas. This force is taken from the Fokker-Planck collision operator 
(cf., Dreicer, 1959) and is given. for a test electron of velocity _v by 
the expression 


where - 


■D 

4ireAcnA 


m r ir- Z 
e e \ f • 


.2 I m f +m e 


f f {v') dy- 


m. 


v - v 


( 5 ) 


m_ 


2 • 


f f {v) is the field particle velocity distribution function, the sum is 
over all field particles, and the integral over all velocity space. For 
the special case where the field particles have Maxwellian distribution 
functions and consist of electrons and protons, equation (5) reduces to. 
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-D 


-e E, 


u 


2A(u) + A(gpU) 


where 


m / m % 

.E - n e/e r a , 3 n ~ ~~ 
c e e p \ rn e 


a e = 




2kT r 


A(u) = <f>(u) - 11 $' (u) 


( 6 ) 


4>(u) = / exp(-t 2 )dt , 

o 

and n and T are the field particle density and temperature respectively. 
In Table II we list the absolute value of F Q , normalized to eE c , as a 
function of u. We see that the dynamical friction increases sharply up 
to the proton thermal speed (6 p u = 1) and then decreases for larger u. 

p 

For high speed test particles (i.e., u » 1) | Fpj decreases as 1/u 
while for those moving at very low speeds (i.e., 3 p u « 1) j F p j becomes 
proportional to u. This behavior of the dynamical friction force is 
responsible for runaway effects associated with high speed particles and 
the drifting of thermal particles in a plasma subject to an externally 
applied field. For a gas with a temperature gradient, we will find that 
we can describe velocity space in terms of a collision-dominated part 
populated by particles which are cooled by dynamical friction as they 
move from higher temperatures to lower temperatures; and a collisionless 
regime populated by high speed test particles which can penetrate down 
through a temperature scale length without cooling substantially. An 
electric field applied along the temperature gradient will have the 
effect of accelerating electrons into the collisionless regime, leading' 
to an increase in the population of the high energy tail. We will 
illustrate these effects more clearly by computing the orbits of high 
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speed electrons in velocity space, with an analysis similar to that used 
by Dreicer (1960). 

Consider a single test electron with velocity v moving through a 
plane-parallel layer composed of electrons and protons whose distribution 
functions are Maxwellian and described by a density, n, and temperature, 
T. For generality we include the effects of an externally applied 
electric field and allow the temperature of the field particles to vary 
along the vertical coordinate, z, of the, layer but, for simplicity, we 
assume that the temperature gradient and density remain constant. Under 
these conditions, the Langevin equation for the change in velocity vis, 
cf„, Chandrasekhar (1943) for details. 


dv 3eE 

dt m a„ 
e e 



-eE 

~ + A(t) 
rn - 
e 


( 7 ) 


In this equation, the acceleration due to particle-particle interactions 
is separated into a time-averaged part, the dynamical friction, and a 
part, A(t), describing fluctuations about the average. The second term 
on the left hand side of equation (7) is simply the high velocity limit 
of the dynamical friction. The parameter, t, represents the time as we 
follow the test electron along its trajectory through the gas, and 
is related to the coordinate, z, by the equation. 


dz _ 
dt 


( 8 ) 


Ignoring the velocity fluctuations produced' by the acceleration A(t) 
(See Dreicer, 1960} and combining equations (7) and (8) we can write 
the Langevin equation as. 
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dT dy t 3e£ c v 
V dz dT a e m e v 3 


-e£ 

m_ 


(9) 


We can reduce equation (9) into two component equations which in non- 
dimensional form are. 


y dx. 

y d0 


u ^y- 
y de 


A 

i 

s 


~2 + P 


n - p 2 ) 


2u, 


( 10 ) 


(ID 


where 


x 2 *uX- , “b = -3Ec /E • 


-1 


e=T/T o . • L = ■ 

and a subscript 2 ero indicates. that the parameter is to be evaluated at 
the initial position of the test electron. The flow lines for. fast 
electrons in velocity space in any plane containing E, are obtained from 
equations (10) and (11) , 


dx 

dp 


o-ift 


1 + px' 


2 


( 12 ) 


The solution to this equation is. 


x 2 „ + „K_ 

x ! + y 1-P 2 


(13) 


where K - ((1+P 0 )x‘ + 2)(1 - p Q ), 


and x o and P Q are the values of x.andp at the initial position of the 
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test electron. Note that equation (13) generates essentially the same 
flow lines obtained by Dreicer. 

The dynamical friction acts to accelerate test particles -in a 
direction which always opposes their motion. This arises simply from 
the fact that the field particles are distributed symmetrically in 
velocity space. For the special case E = 0 the flow lines are generated 
by the equation p = constant so that, in velocity space, the test particles 
simply follow straight line orbits leading into the origin. The relevant 
problem for a test particle moving through a finite layer is to find its 
final velocity after it leaves the layer, given its initial velocity u 0 . 

By computing the work done by the dynamical friction on a test electron 
as it moves through the layer and given its initial velocity, we can 
derive its final velocity after it leaves the layer. As shown in Appendix 1, 
this is given in general by the equation, 



where u^ and are the final and initial speeds normalized to the 
to / ^ \ 


thermal speed at the top of the layer; S c = 




is the work done 


(normalized to kT o ) by the dynamical friction on a thermal electron 
(at T ), which travels a temperature scale length (L) along the temperature 
gradient; 0 Q andO^ are the initial and final temperatures normalized 
to T . Equation (14) neglects the effect of deflections in slowing a 
beam of test particles. 


In Figure 7, we plot, on a v vs. v graph the final velocities for 
test electrons injected into a given layer in the transition region 
(characterized by its critical speed, u c , computed from Dupree's model) 
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2 2 2 2 * . . 

for two initial energies, u = 1.5u and u = 2.0u . The circle in 

U v U U 

The critical speed (u c ) referred to here is that for a 90° deflection 


this plot represents the distribution of initial velocities, all with 
the same magnitude, for test electrons injected at different angles to 
the temperature gradient. Note that those particles with positive 
initial y-values are injected at the bottom of the layer; while those 
with negative u- values are injected at the top. The straight lines are 
the particle flow lines which are drawn in to simplify the tracing of 
the particle trajectory from its initial velocity to its final velocity • 
after it leaves the layer. Clearly, particles injected at large angles 
(small y) to the temperature gradient never make it out of the layer. 

•k 

These test electrons simply come to rest somewhere in the layer itself. 


*We note that a single test particle can never actually come to rest 
because of the effect of the velocity fluctuations, A(t), which we ignored 
in this analysis. Furthermore, remember that our analysis only applies as 
long as the test particle speed is larger than the local thermal speed. 


On the other hand, those electrons moving along the temperature gradient 
have small changes in their initial speeds. Another interesting effect 
is associated with the y dependence of the final speeds. We see that 
the relative change in particle energies does not vary substantially 
with y until a critical value is reached, beyond which the test particle 
decelerates rapidly. For higher initial energies, this critical y-value 
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decreases. Thus, even if a test particle undergoes fluctuations in 
velocity space about its mean orbit, it does- not slow substantially 
until it reaches large angles to the temperature gradient. 

The effect of the statistical fluctuations in velocity, A(t), which 
v/as not included in this analysis, is to smear out the particle orbits 
in velocity space. In other words, a single test electron can jump from 
one flow line to another in the course of its motion through a given 
layer. The degree to which a particle is deflected to and from various 
orbits depends on its speed. The larger the particle speed, the smaller 
the fluctuations. Clearly, those electrons moving along the temperature 
gradient have the greatest chance of escaping a given layer without a 
substantial change in their initial speeds or directions. Those traveling 
at large angles to the temperature gradient, will slow because of their 
longer path lengths and will also deflect more easily to different flow 
lines. Some will be deflected into a flow line directed along the 
temperature gradient; however, if their speeds have decreased substantially 
by this time, they can be easily deflected back out again before escaping 
the layer. These particles become trapped in the thermal pool. Thus, 
the dynamical friction and A(t) combine to thermalize and isotropize the 
distribution of test particles to the field particle distribution. 

However, they become ineffective for particles at high speeds and 
traveling at small angles with respect to the temperature gradient. 

This is what leads to anisotropies and a high energy tail in velocity 
space for a plasma with a temperature gradient. We can conclude that 
velocity space can be broken up into a thermal part and a collisionless 
part populated by high energy electrons which stream through a temperature 
scale length without altering their energies substantially. 
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A similar analysis can be carried out for the case where an electric 
field is applied along the temperature gradient- Combining equations 
(11) and (13), we obtain a transcendental equation for the final p-value 
of the test electron; namely. 



2 

1+P 





1 + 


2 2 
Vo 



(15) 


From this equation and the particle trajectories (equation 13), we 

obtain the final test particle velocities. Equation (15) was solved 

numerically with an iteration scheme. Our results are displayed in 

Figures 8-11 with a format similar to that used for zero electric 

field. In this case, however, the results depend on both the value of 

the critical speed, u c , for a given layer in the transition region and 

on 3 = E/E c - In Figure 8 , the final velocities are associated with 

2 2 2 2 2 

initial energies of u = 1.5u and u = 2.0u where u is the minimum 

U U w v 

, o 

value of (u ) D plotted in Figure 4 and 3 was taken to be 0.5. In Figure 9, 

the initial velocities are the same as in Figure 8 , however,' the value 

of 3 was changed to 1.0-. We see that the larger the initial velocities 

the smaller the final relative change in velocities. In addition, if 

the electric field is increased the relative change in velocities increases, 

2 

In Figures 10-11, we plot the final velocities for a value of u c equal 

2 

to twice the minimum value of (u c )p for the transition region and for 
the same values of 3. Comparing these two sets of figures, we find that 
the acceleration due to the electric field increases relative to that 
due to the dynamical friction as we move to higher speeds. This arises 
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2 

from the fact that the dynamical friction force falls off as 1/u while 
the electric field is independent of velocity. 

The most interesting features of these plots are the particle 
orbits. As pointed out by Dreicer (1960), there exists a critical 
surface which- separates velocity space into a region where dynamical 
friction dominates and a part where the electric field dominates. The 
minimum distance of this surface from the origin occurs at ^ = - 1 and 
is given as u = u^ where u^ = 3E C /E. The larger the electric field, the 
smaller the particle speed needed for the electric field to dominate 
over dynamical friction. Thus the region in velocity space which is 
effectively co.llisionless becomes larger with increasing E. We also 
observe that this critical speed increases as y decreases. What is more 
important, however, is the behavior of the final velocities as a function 
of y. We find that the relative change in velocities remains approximately 
constant up to a critical angle, below . which a test particle experiences 
larger changes in velocity. Indeed, the latter region is quite narrow 
in y and becomes narrower the larger the initial speed of the test 
particles. 

Finally, the effect of an electric field is essentially to accelerate 
test particles into the collisionless regime. If the electric field 
points along the temperature gradient, then the asymmetry produced by 
the temperature gradient is enhanced i.e.,'the number of high energy 
electrons moving down the temperature gradient is increased. The opposite 
applies if the electric field points in a direction opposite to the 
temperature gradient. Furthermore, the electric field introduces a 
critical surface within which a particle's motion is dominated by 
dynamical friction while outside of which the motion is dominated by the 


electric field. 
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4. ESTIMATE OF THE ELECTRON DISTRIBUTION FUNCTION 
We have shown that the electron distribution function for velocities 
parallel or anti -parallel to a temperature gradient could be divided 
into a collision-dominated regime and a collisionless regime. If the 
collision dominated part of the distribution function can be represented 
by a Maxwellian, then for all v x , and for a steady state with a 
temperature gradient in the z-direction. 


f e ( vyv z) = n o 


m„ 


2vkT. 


3/2 


exp 


m„ 


2kT, 


(v 2 +v 2 +v 2 ) 
v x y z‘ 


(16) 


for v z >_ -v zc v/hile for v z £ v zc 


m_ 


f e (v x’V V z) = ST £X P 


m. 


2kT 


(v 2 +v 2 ) 
' x y' 


f tail ( V z) 


where v is the 'critical 1 speed beyond which the electrons become 

collisionless. We note from equation (16) that the tail is attached 

only at negative z-velocities. This stems from the argument that there 

are very few collisionless test particles originating from the lower 

temperatures and very few high energy particles, from high temperatures, 

deflected back up the temperature gradient. 

Our objective, then, is to estimate the velocity dependence of the 

tail (f tal --j [v z ,z]) of the electron distribution function throughout the 

transition region, given the density and temperature profiles from 

Dupree's model. We note first that the tail at a temperature, T, must 

ultimately originate from the thermal part of distribution functions at 

* 

higher temperatures. As a result, the contribution to the tail at 


This conclusion stems primarily from the fact that the corona is isothermal 
and that the tail is populated by collisionless electrons which stream to 


lower energies without changing their initial energies.- 
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temperature, T, and velocity, v 2 , from the distribution function at 
T + AT and velocity v z is given by. 


f tan< v r T > = "( T+AT > 

[, Vz 1 
exp [“ 2irk(T+AT ) J “ v zc 




27rk(T+AT) / 


h 


v, r <T+AT)‘ < v z < -v zc (T) (‘1 7 ) 


where the upper limit on v z arises .from the fact that particles with 

larger velocities are collision dominated (not part of the tail) and the 

lower limit is a result of the fact that particles at lower velocities 

are also collisionless-at T + AT and therefore originated ultimately 

from the thermal part of the distribution functions with temperatures 

greater than T + AT. If we take the limit of equation (17) as AT 

approaches zero, we see that each temperature, T, contributes, at a 

single velocity v z = v (T), an amount given by a Maxwellian evaluated 

at v z » to the tail at lower temperatures. Therefore, the tail at temperature, 

T, is given .in general by. 


Wv T > ■ n < r > 


1 \ % 


2 , 

m e . 

exp 

3 

< 

N 

1 2irkT‘ i. 

L 2kT ‘ J 


08 ) 


for v z £ -v zc (T) and where T' is a function of v , derived from the 
equation. 


v zc {r) = v z ’ (19) . 

The form for the tail is complete, given the critical speed as a 
function of temperature. In Section 2, we estimated the critical energies 
beyond which electrons become collisionless from ratios of appropriate 



- 20 - 


time scales and from collision-dominated theory. Three sets of critical 

2 2 

energies were obtained, (u c ) D {90° deflection), (u c ) £ (thermal izati on) 

o 

and u (collision-dominated theory). Since the latter is based on a 
c 

precise mathematical solution of the Boltzmann equation, we feel that 

u 2 gives a good indication for the location of the tail in velocity space. 
2 

However, u was also derived assuming that the first order distribution 
c 

function f^ was equal in magnitude to f^. Since strong departures 

from Maxwellian will occur even if f^ is one-tenth of f^, the critical 

2 

energy should be chosen to be less than u . We have chosen to work with 
? 2 

(u )r: which is always less than u^ throughout the transition region, 
c t c 

This choice is somewhat arbitrary and should be considered a lower limit 

2 

on the critical energy since electrons with energy less than (u c )^ are 
thermal i zed. We will also present results for a critical energy equal to 

JL. O 

(0.4r(u£)p. This particular choice arises in connection with the calcu- 
lations of particle orbits presented in Section 3 (see also Roussel-Dupre, 
1979 5 p. 97) and should be' considered an upper limit to the critical 
energy. I 

2 

Given the results of Section 2 for the critical energies (u^) F , we 


have 



= 3.28 x. 10 2 


n£nA 

f dT' 1 

-1 

T 1 

. 

[dz j 



1/3 


( 20 ) 


Combining equations (18), (19), and (20) and assuming that the gas 
pressure and the quantity q‘ = T ,5//2 dT'/dh are constant, we obtain 


f tail<V T o> = n c 


nr 

e 

2irkT 


k 


I .2 


9/7 


zc 


x e-xp 


-u 2/7 u 12/7 

z zc 


( 21 ) 


for u z < -u zc . 
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Equation (21) gives the final form for the tail of the distribution 
function through the transition region. In Figure 12, we plot the ratio 
of f + ,--i(u v ,T„) to a Maxwellian evaluated at u. as a function of u_/u__ 

la I I Z 0 Z Z ZC 

2 

for several values of u . We find that this ratio increases rapidly 
with velocity i.e., the tail is overpopulated compared to a Maxwellian. 

Finally, implicit to the derivation of the velocity dependence" of 
the tail was the assumption that the collision! ess particles gain or 
lose only a small fraction of their initial energy over the extent of 
the transition region as a result of work done on them by external 
forces such as gravity. This is an excellent approximation since the 
transition region is so thin and since these collisionless particles 
have very high energies. In addition, the assumption that q' is constant 
is also a very good approximation for the transition region since most, 
if not all of the collisionless particles originate from the region 
T = 10 5 * 2 - 10 6 * 2 K where T 5y/2 dT/dh is indeed a constant. 

5. MOMENTS OF THE ELECTRON DISTRIBUTION FUNCTION 

In deriving equation (18) we assumed that the thermal part of the 
distribution function could be approximated by a Maxwellian. In fact, 
the actual distribution function is the sum of a Maxwellian plus small 
correction terms which result from nonuniformities in the density, 
temperature and mean flow of the gas. If the correction terms are small 
they may be computed from a perturbation analysis - the Chapman-Enskog theory. 
We found in Section 2 that the first order correction term computed by 
Spitzer and Harm (1953) for electrons in a nonuniform gas, is small for 
electron energies less than a critical energy, u c , plotted as a function 
of temperature in the transition region in Figure 2. Since the correction 
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term is small compared to a Maxwellian for speeds defined to be in the 

thermal part of the distribution function, we were justified in our 

approximation of the previous section. However, we cannot completely 

ignore the contribution of the first order term to the overall distribution 

function since its nonuniform nature in velocity space contributes to 

the moments which determine the mean gas flow and heat flux at a given 

height in the transition region. Since perturbation theory breaks down 

beyond the critical speed, u c , this correction term can only apply to 

the collision dominated part of the distribution function and, therefore, 

does not affect the tail. In fact, the tail was computed because collision 

dominated theory broke down at high velocities. The total electron 

distribution function for all v and v becomes, 

x y 


f e<vyv z) = f e ! 0 )( WV z) + f e < 1 )( VW z) 


for v ~ v 7r and 


( 22 ) 


( 0 ) 


f e ( Vy v z’ z > = f e (v x’V z) f tail ( V z) 


for v _< -v where f ^ is a Maxwellian; f ^ is taken from Spitzer 
and Harm (1953) and is given by equation (21). 

The electron density, mean flow temperature, and heat flux in the 
transition region are defined in terms of the electron distribution 
function as: 


n e ( z ) = J 

v e (z) = J 

m 

Temperature: T(z) = 


Density: 
Mean Flow: 


dv f e (v,z) 
dv v f e (v,z) 


: e ) f e {v -' z) 


J dv (y-v £ 

q e (z) = J dv (v-v e ) 2 v f 0 (v,z) 


Heat Flux: 
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From equation (19), we then find, 


"e' 2 ' * "o 


.12 


V 2 (1 + *(„ c )) + Z£- festgL 


dt 


(23) 


= 2 e< z > ‘ 


4 , , exp(-u c ) 

J7¥ + 2 A" 


x (7 ,exp(u c 2 ) u 2 E 3 (u|) - 1) 


(24) 


V z >=i; T „l 4 - 3 * < u c» 



1 2kT 1 

% " 

q (z) = -n kT 
M ze oo 

0 


i m e 1 



A w 


3/ir 
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7 exp(-u c ) 
/it 


In.Z 


3u 


c , 27 . 6 . 12 . 12 

- 7 * 7 . 

c c c / 


(26) 


where <j> is the error function; is the third exponential integral; 
u is the critical speed at 'which the tail is attached; and u__ is the 

v* * Z 6 

mean flow of the electron gas (which has a z-component only) normalized 
to the local temperature, T , and I (x) is defined by. 


1 » W 'I 


t n exp(-t 2 ) D g (t) dt. 


We note that n g and T e reduce to n Q and T q as u c approaches infinity. 

This must be the case since the distribution function then reduces to. a 

Maxwellian plus a correction term which does not contribute to either of 

these moments (f ^ « p). As u„ approaches infinity the mean flow and 
6 c 

heat flux also. reduce to their values predicted from collision dominated 
theory. 
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We can evaluate the effect which the tail of the distribution 
function has on the various moments discussed above by comparing the 
actual moments to the values predicted by collision dominated theory 

n -n T -T 

(i.e., u c -»- ~). In Table 3, we list values for — , — and 




as a function of temperature in the transition region. We find 


that the tail has very little effect on the temperature and density of 
the gas as predicted from collision-dominated theory. In other words, 
n o and T q represent the total density and mean kinetic energy of the gas 
to a high degree of accuracy. On the other hand, the heat flux is 
enhanced from that predicted by collision dominated theory by as much as 
thirty percent where the temperature gradient is steepest. 

In the case of the mean flow for the electrons, we cannot define a 
u u 

quantity— 6 —— since, assuming that the ions and protons are 
' u 

zo 

Effectively stationary relative to the electrons and that no net charge 
enters or leaves the transition region, the electric current must equal 
zero or u ze and u zp equal zero. In the presence of a temperature 
gradient, zero net current is maintained by -a self-consistent electric 
field which drives cold electrons up the temperature gradient to compensate 
for the flow of hot electrons down the temperature gradient. We can 
compute the electric field, E, needed to maintain zero net current by 
setting equation (24) to zero. In the collison dominated case (u -*-«>}, 
where the mean ionic charge (Z) is equal to one, this electric field is 
given by, 

= -0.703 £ S- 
o e dz 


In Table 3, we list values for (E - E Q )/E o as a function of temperature 
in the transition region. We see that the tail does not affect the 
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value of the collision dominated electric field to any significant 
degree. We also list in Table 4, values for E/E as a function of 

V 

temperature; we find that the self-consistent electric field is small 
compared to the dynamical friction. Correspondingly, it is still a good 
approximation to assume that the high energy particles streaming down 
the temperature gradient do not alter their energies substantially 
because of acceleration or deceleration due to an external force (See 
the results of Section 4). 

Because of the difficulty in determining a precise value for the 
critical energy we have presented calculations for minimum and maximum 
values. The results presented in Table 4 correspond to a minimum value 
for the location of the tail in velocity space. In Table 4, we present 


values for 




e ^o 


computed with the maximum critical energy given by 


(0.4)^ (u ) D . We see that these values are all negative, meaning that 

collision-dominated theory predicts a larger heat flux than the distribution 

function with a tail; however, this is not physically valid since particles 

cannot transport energy any faster than if they stream freely as in the 

case of the high energy tail. These results indicate that if the critical 

energies (used to obtain the results in Table 4) at which the tail is 

attached are correct then collision-dominated theory tends to overestimate 

the electron heat flux. This is not surprising since as shown in 

Section 2, collision-dominated theory breaks down beyond a critical 

energy plotted for the transition region in Figure 2. On the other 

hand, the results obtained with the critical energy for the tail given 

p p 

by u c = (u c ) E , give heat fluxes which are greater than predicted by 
collision-dominated theory. We can conclude that the actual heat flux 
is quite sensitive to the location of the tail in velocity space. These 
results also suggest that the critical energy for the location of the 
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tail cannot exceed the critical energy beyond which collision-dominated 

* 

theory breaks down. We note that the latter energies are closer to the 

This suggestion stems from the fact that we have incorporated the first 
order distribution function, f^ s in our total distribution function. 

minimum critical energies used in our calculations. Indeed, as already 
suggested, the minimum critical energies, (u c ) E , are probably a better 
estimate of the location of the high energy tails. An exact solution 
to this problem would require solving the Boltzmann-Fokker-Planck equation 
numerically - a task of enormous proportions. 

6. SUMMARY 

In Section 2 of this paper, we found that, based on a comparison of 
appropriate time scales, the electron and proton distribution functions 
were non-Maxwell i an in the solar transition region in the sense that a 
high energy tail composed of hot electrons and protons streaming down 
from the corona would result. We also found that this anisotropy in the 
proton and electron distribution functions would not be reflected by the 
ions since their higher charge causes them to be collision-dominated and 
results in nearly Maxwellian distribution functions. In addition, we 
pointed out that the proton and electron temperatures need not be equal 
through the transition region and that the ions would reflect the 
proton temperature. If the proton temperature differs from the electron 
temperature in the corona then this condition would persist through the 
transition region. In a steady state, however, the latter condition 
would necessitate a source or sink of energy for one or the other 
species in the corona. We suggest that this would represent an interesting 
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line for future research since observations show line widths which are 
much larger than the thermal widths expected from electron temperatures 
derived from ionization equilibrium calculations. 

In Section 3, we computed velocity space orbits for test electrons 
moving through a layer composed of field particles with Maxwellian 
distribution functions characterized by a temperature which varies 
through the layer. In addition, we computed the final speeds for these 
test electrons after they leave a layer whose depth is defined by a 
temperature scale height. The results were obtained from the Langevin 
equation which includes the effects of a dynamical friction force and an 
externally applied field. For the case of zero electric field, the 
velocity space orbits are straight lines (given by p = constant) leading 
into the origin in velocity space. In computing the final energies for 
the test electrons, we found that an electron's energy does not change 
substantially over a temperature scale length for energies greater than 
a critical energy and that this result applies over a broad range of 
y-values up to a critical value which is a function of the particle's 
energy. When an external electric field was included in our calculations 
we found its effect was essentially to accelerate test electrons into 
the collisionless regime. If the electric field points along the tempera 
ture gradient, then the asymmetry produced by the temperature gradient 
is enhanced i.e. , the number of high energy electrons moving down the 
temperature gradient is increased. The opposite applies if the electric 
field points down the temperature gradient. Furthermore, as shown by 
Dreicer (I960),, the electric field introduces a critical surface in 
velocity space. If a particle's energy is less than the energy along 
this surface, then the dynamical friction dominates the particle's 
motion. If it is greater, then the particle's motion is dominated by 
the. electric field. 
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The main conclusion of Section 3 was that, for a plasma with a 
temperature gradient, velocity space can be broken up into a thermal 
part, and a collisionless part populated by high energy electrons which 
stream through a temperature scale length without altering their energies 
substantially. The net result is that the electron distribution functions 
are characterized by high energy tails. In Section 4 we estimated 
the magnitude of these tails, for conditions appropriate to the transition 
region, and found that they were strongly over-populated relative to a 
Maxwellian. 

In Section 5, we used the results of Section 4 to recompute various 
moments of the electron distribution functions. We found that the tail 
had a negligible effect on the total local density of electrons, the 
electron temperature, or on the critical self-consistent electric field, 
needed to maintain zero net current. On the other hand, the effect of 
the high energy tail is to enhance the electron heat flux over that com- 
puted from collision-dominated theory by a maximum of 30% in the temperature 
c 2 S 4 

range from 10 ’ - 10 K. We also found, however, that these results 

are sensitive to the location of the high energy tail in velocity space 
and that detailed calculations, which would involve solving the Boltzmann- 
Fokker-Planck equation numerically, are necessary in order to determine 
the exact magnitude of these effects. Nevertheless, we feel that the 
minimum critical energies represent a good approximation for the locations 
of the tails and that the results obtained with these critical energies 
are reasonable approximations for the enhancement in heat flux. 

The enhancement in the heat flux, which we computed for Dupree's 
model will not have a serious effect on the energy budget of the corona 
and the transition region. However, it is important to realize that the 
computed enhancement is sensitive to the location of the critical energy. 
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Indeed, the increase in heat flux relative to that computed from collision- 
dominated theory may be significant for regions which have lower densities 
and/or larger temperature gradients (i.e., lower critical energies} than 
obtained from Dupree's model. A similar analysis for coronal holes and 
active regions, for example, may yield interesting results.' 
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Table I 

2 

Fraction f of Particles with Energy Above u 
in a Maxwellian Distribution 


u 2 

1 - f 

f 

1.00(-4) 

7,S2{-7) 

1.00 

1.00C-2) 

7.48(-4) 

1.00 

l.OO(-l) 

2-32 (-2) 

9.77 (-1) 

5.00(-l) 

2.01 (-1) 

7.99(-l ) 

1.00 

4.28(-l ) 

5.72(-l) 

2.00 

7.36(— 1 ) 

2.64(-l) 

4.00 

9.54(-l ) 

4. 60( -2) 

6.00 

9.93(-l ) 

6.85(-3) 

8.00 

9.99(-l) 

1.07 (-3) 

1.00(1) 

1.00 

1 - 62( -4 ) 

5.00(1) 

1.00 

1.54(-21) 

1.00(2) 

1.00 

4.20(-43) 




i t 
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Table II 

The Dynamical Friction Force for Electrons 


u 

F D / e E c 

0.0 . 

0.00 

1.0 (- 4 ) 

5.90 

1.0 (- 2 ) 

5.36 ( 2 ) 

2.0 (- 2 ) 

7.83 ( 2 ) 

2.3 (- 2 ) 

7.84 ( 2 ) 

3.0 (- 2 ) 

7.37 ( 2 ) 

1.0 (- 1 ) 

1.00 ( 2 ) 

2.0 (- 1 ) 

2.53 ( 1 ) 

3.0 

1.15 ( 1 ) 

5.0 (- 1 ) 

4.65 

1.0 

1.86 

1.4 

1.26 

1.8 

8.70 (- 1 ) 

2.0 

7.27 (- 1 ) 

2.5 

4.80 (- 1 ) 

3.0 

3.33 (- 1 ) 

4.0 

1.88 (- 1 ) 

5.0 

1.20 (- 1 ) 
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Table III 



Moments of the Electron 

Distribution 

Functions 



Log T 

n - n Q 
n o 

T - T o 

T o 

% ’ q o 
% 

E - E o 
E o 

E 

4.7 

7. 56 (-11 } 

4.4(-9) 

-1.00(-5) 

2 . 03 (—6) 

1.31 (-4) 

4.8 

T.88(-6) 

7.2(-5). 

-5.11 (-3) 

-1.52(-3) 

‘ 7.89(-4) 

5.0 

4.09(-4) 

1.2(-2) 

1 . 21 ( -1 ) 

1.20(-2) 

4.82(-3) 

5.2 

1.94 (-3) 

5.1 (-2) 

2.83(-l) 

3.50(-3) 

1.12(-2) 

5.4 

1.31 (-3) 

3.6(-2) 

2.92(-l) 

1.38(-2) 

8.87 (-3 ) 

5.6 

7. 37 (-4) 

2.1 (-2) 

1 .69{-l ) 

1.44{-2) 

6. 52 (-3) 

5.8 

.4.09(-4) 

1.2(-2) 

8.59(-2) 

5.98(-4) 

5. 02 (-3) 

6.0 

2.05(-4) 

6.2(— 3) 

4.73'(-2) 

-2.29(-3) 

3.66(-3) 

6.2 

1 .13(-4) 

3 . 5(— 3 ) 

1 -89(— 2) 

3 . 32( — 3 ) 

2.87 (—3) 


The subscript (o) refers to moments computed from collision 
dominated theory. Those without the subscript are computed 
with the effect of the tail included. 
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Table IV 

St 

Comparison of Electron Heat Flux Calculations 


Log T 

2 

U 2C 

% 

4.7 

66 

0 

4.8 

27 

-7.70(-7) 

5.0 

11 

-3.45(-2) 

5.2 

7.0 

-1.80(-1) 

5.4 

8.2 

-1.02(-1) 

5.6 

9.5 

-7.04(-2) 

5.8 

11 

-3.51 (-2) 

6.0 

13 

-1.41 (-2) 

6.2 

14 

-5.45(-3) 


* 


is the heat flux computed with the effect of the of the tail 
deluded while q Q is the collision dominated value. 



FIGURE CAPTIONS 


Figure 1 Ratio of electron mean free path to temperature scale length 
for thermal electrons in the Solar Transition Region. 

Figure 2 Critical energies at which the first order electron distribution 
function equals the Maxwellian. 

Figure 3 Critical energies at which the first order proton distribution - 
function equals the Maxwellian. 

Figure 4 Critical energies for 90° deflection and thermal izati on of 
electrons. 

Figure 5 Critical energies for a 90° deflection and theramal izati on of 
protons. 

Figure 6 Critical energies for a 90° deflection and thermal izati on of 
Si IV. 

Figure 7 Final velocities for test electrons injected into a given 
layer in the transition region for two initial energies. 

Figures 8-11 Final velocities for test electrons injected into a 
given layer in the transition region with an electric field. 

Figure 12 Ratio of the electron high energy tail to a Maxwellian. 
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APPENDIX I 

DERIVATION OF THE WORK DONE ON A FAST ELECTRON OVER 
A TEMPERATURE SCALE LENGTH 


The dynamical friction force which appears in the Langevin equation 
includes the effect of deflections in slowing a beam of particles; The 
deflections, however, do not remove energy from the beam but singly 
redirect some of that energy into the perpendicular direction. For our 
purposes, we are interested in computing the amount of energy lost by a 
fast. electron to the field particles during Coulomb encounters over a 
temperature scale length. More precisely, given a test electron's 
initial energy, we wish to compute its final energy after it travels a 
temperature scale length. 

The change in energy experienced by a test particle during a single 
encounter can be written as. 


AE = j (Av^ + A \} jJ- + 2v ^ //) 


(i-n 


where perpendicular and parallel here are taken relative to the initial 
direction of motion. Averaging (Irl) over allimpact parameters and 
field particle velocities, assuming a Maxwellian distribution function 
for the field particles, we obtain a general expression for the average 
rate of energy exchange with the field particles (See Spitzer 1962), 


<AE> 

kT 


2n t r t a t 3/^ 1 ( * (6fU) 


n t Z t 


- (i + e f - 2 )A (B fU )) 


(I-?) 
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In the absence of any external forces, the rate of energy loss for a 
fast electron moving through an electron-proton gas is obtained from 
equation (1-2) and is given by. 


d 



dt 


2n r a a 
e e e 


(1-3) 


where we have assumed that n g = n . The time, t, in the above equation 
is taken along a particle trajectory and is related to the 'z-coordinate 
in our layer by. 


dz _ 
dt 


W. 


(1-4) 


We can now make a change of variables where the temperature along a 
particle's trajectory replaces the time in equation (1-3). Assuming 
that the field particle densities and the temperature gradient are 
constant, we have, 


d(u 2 e) 2 

de 


T / p 

where 6 = ; T and a„ = 

o o 2kT 


The solution to (1-5) is simply, 




(1-5) 


( 1 - 6 ) 





A-3 


where is the final speed normalized to the thermal speed at 
T q and u Q is the initial particle speeds. Given that, 

2/1 d"r\ 

^ n e r e a o (~ dz" ) = e( l lja ' t ' ion U - 6) reduces to the desired 

result - equation (14), 



